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ABSTRACT 

3 : 

Here we discuss the X-ray emission properties from the hot thermalized 
^ ■ plasma that results from the collisions of individual stellar winds and supernovae 

ejecta within rich and compact star clusters. We propose a simple analytical way 
of estimating the X-ray emission generated by super star clusters and derive an 
expression that indicates how this X-ray emission depends on the main cluster 
. parameters. Our model predicts that the X-ray luminosity from the star cluster 

region is highly dependent on the star cluster wind terminal speed, a quantity 
related to the temperature of the thermalized ejecta. We have also compared 
the X-ray luminosity from the SSC plasma with the luminosity of the interstel- 
lar bubbles generated from the mechanical interaction of the high velocity star 
q . cluster winds with the ISM. We found that the hard (2.0 keV - 8.0 keV) X-ray 

£j . emission is usually dominated by the hotter SSC plasma whereas the soft (0.3 

keV - 2.0 keV) component is dominated by the bubble plasma. This implies that 
£> . compact and massive star clusters should be detected as point-like hard X-ray 

•rH . 

^ . sources embedded into extended regions of soft diffuse X-ray emission. We also 

compared our results with predictions from the population synthesis models that 
take into consideration binary systems and found that in the case of young, mas- 
sive and compact super star clusters the X-ray emission from the thermalized 
star cluster plasma may be comparable or even larger than that expected from 
the HMXB population. 

Subject headings: Star clusters: winds, bubbles, X-ray emission 



1. Introduction 



In many starburst and in interacting and merging galaxies a substantial fraction of the 
star formation is concentrated in a number of compact, young and massive stellar clusters 
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or super stellar clusters (SSCs; see, for example, Holtzman et al. 1992; Ho, 1997; O'Connell, 
2004 and references therein). The high stellar densities, the large energy and mass deposition 
rates, provided by stellar winds and supernovae explosions (SNe), suggest that SSCs are 
potentially strong X-ray emitters (Chevalier, 1992). 

Indeed the X-ray emission from NGC 3603, Arches cluster in our Galaxy and R136 
in the LMC, the local analogies to low mass SSCs, has been detected (Moffat et al, 2002; 
Yusef-Zadeh et al. 2002; Stevens & Hartwell 2003). The X-ray emission from Arches cluster 
has a diffuse component and five discrete sources. Overall, 40% of the observed emission is 
diffuse and 60% comes from the discrete sources. The absorption-corrected X-ray luminosity 
(between 0.5 - 8 keV) is Lx ~ (0.5 — 3) x 10 35 erg s -1 and the best spectral fit requires the 
metallicity of the X-ray plasma to be 4 - 5 times solar (Law & Yusef-Zadeh, 2004). Similarly, 
the X-ray luminosity from R136 in 30 Doradus is L x = 5.5 x 10 34 erg s _1 (Stevens & Hartwell 
2003). On the other hand, an examination of the XMM- Newton observations led Smith & 
Wang (2004) to the discovery of the diffuse X-ray emission from a 100 pc superbubble, 30 
Doradus C, associated with a rich OB-association LH 90, located south-west of the center 
of 30 Doradus. The examination of the X-ray spectrum indicates also an overabundance of 
a-elements in the bubble plasma by a factor ~ 3, compared to the metallicity of the LMC. 

X-ray emission detected in distant galaxies presents overall two components: an ex- 
tended diffuse component that seems to be related to a collection of superbubbles associated 
to individual star clusters, and a number of unresolved point-like sources (Summers et al. 
2004; Metz et al. 2004; Smith et al. 2005). In some cases, the X-ray emission shows good 
coincidence with the young stellar clusters (Metz et al. 2004). However, there are also many 
examples of X-ray regions that do not present any evidence of active star formation. In 
other cases the X-ray sources seem to be shifted relative to the nearby star clusters and 
some massive clusters seem to emit no X-rays at all (Kaaret et al. 2004). Thus the origin 
of the X-ray emission in distant galaxies and the dependence of the detected X-ray emission 
on the parameters of the embedded star clusters remain controversial. 

Here we discuss how the mechanical energy deposited by stellar clusters restructures the 
host galaxy ISM while producing a hot plasma which should be detected in the X-ray regime 
(see section 2). In sections 3 and 4 we demonstrate that the mechanical energy from SSCs 
leads to a two-component plasma model: a high temperature, overabundant in a— elements 
plasma, ejected from the star cluster, and a lower temperature, extended, interstellar bubble 
component. Here we discuss how the luminosity of the hot component depends on the star 
cluster parameters and in section 5 we compare the contributions from both components 
when evolving in different interstellar environments and as a function of time. 
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2. Super Star Cluster Winds and their Impact on the ISM 

The interaction of the high velocity outflow that results from the energy injected by 
multiple stellar winds and supernovae (SNe) explosions within a star cluster volume, leads to 
a major re-structuring of the ISM. The evolution under the assumption of a uniform density 
case, causes, as in the case of a stellar wind (Weaver et al. 1977), several physically distinct 
regions (see Figure 1), promoted by the shock waves (inner and outer shocks) that result 
from the wind-ISM interaction. The X-ray emission from such structure, must take first 
into account the central star cluster (zone A), filled with the thermalized, high temperature 
plasma, injected by winds and SNe. There the large central overpressure that results from 
thermalization, accelerates the injected gas and eventually blows it out of the star cluster 
volume while generating a stationary wind (Chevalier & Clegg, 1985). The central zone A, 
a powerful X-ray emitter, is surrounded by the free- wind region (zone B in Figure 1) where 
the high metallicity matter emanating from the SSC volume rapidly acquires its terminal 
velocity (iu,e»)j while its density and temperature begin to steadily drop as r~ 2 and r _4//3 . 
Zone B is also, in the case of SSCs, a source of X-ray emission. Its main contribution, 
as shown below, arises from regions in the immediate neighborhood of the SSC, although 
formally X-rays are emitted throughout the whole zone down to its outer boundary: the 
reverse inner shock. An exception to this latter issue is for massive and compact SSCs whose 
winds become strongly radiative, bringing the temperature well below the lower limit for 
X-ray emission (T cut = 5 x 10 5 K), long before reaching the inner shock (Silich et al. 2003, 
2004). Whether the SSC wind undergoes cooling or not, it eventually crosses the inner shock 
to be once again thermalized. This generates and maintains the hot superbubble (zone C), 
the reservoir of thermal energy that drives the outer shock into the ISM. The hot superbubble 
spans from the location of the inner shock to the contact discontinuity that separates the 
hot wind from the expanding outer shell of interstellar swept up matter, accelerated by 
the leading shock (zone D). Given the strength of the shocks, at first the two outer zones 
are strong X-ray emitters. Soon however, particularly for large ISM densitites (tiism > 1 
cm -3 ), radiative cooling brings the temperature down in zone D, causing its condensation 
into a thin, dense and cool expanding outer shell. In cases in which magnetic fields are not 
considered the internal structure of zone C, its chemical composition and its X-ray emissivity, 
depend strongly on the evolution time, t, and on the rate at which mass is evaporated from 
the cold outer shell into the hot superbubble interior (see Bisnovatyi-Kogan & Silich, 1995 
and references therein). 
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Fig. 1. — SSC winds and the structure of the ISM. Schematic representation of the structure 
generated by the mass and energy deposition rates within SSCs. The central zone (A) 
represents the SSC volume itself, where stellar winds and SNe from the evolving cluster are 
thoroughly thermalized, causing the large thermal pressure that drives the SSC wind. The 
remaining concentric zones are, as in Weaver et al. 1977, (B) the free-wind region, (C) the 
hot superbubble and (D) the shell of swept up matter, all surrounded by the unperturbed 
ISM. Note that zones A and B produce an insignificant X-ray emission in the case of winds 
produced by single stars (see Weaver et al 1977). 
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3. The Inner Structure of the X-ray region 

In the adiabatic solution the temperature and density of the gas at the star cluster 
surface, r = Rsc, may be inferred from the mass and energy conservation laws (see equations 
4 and 6 in Silich et al. 2004). They are: 

_ 2( 7 - 1)^ L sc 7-1 »i V 2 11v1 nV 2 k m 

MSC (l + 1\ ^SC ._ o o „ m-24 L 38 „ __ 3 



±<kR 2 sc c S c \7-l/ 27ri?| c Vi i0O ^lc,i y iooo 

where L^c* and M^c* are the mechanical energy and mass deposition rates, k is the Boltzmann 
constant, 7 = 5/3 is the ratio of specific heats, /ij = 14/23m# is the mean mass per particle 
for a fully ionized plasma, c S c = {ikT '/ fii) 1 ^ 2 = [(7 — l)/(7+l)] 1 / 2 Va i0O , is the sound velocity 
at the star cluster edge, and Va,oo = (ZLsc/Msc) 1 ^ 2 , is the resultant wind adiabatic terminal 
speed, L 38 is the mechanical luminosity of the SSC measured in units of 10 38 erg s^ 1 , Rsc,i 
is the star cluster radius measured in units of 1 pc, and V1000 is the wind adiabatic terminal 
speed in units of 1000 km s -1 . The temperature and density of the gas inside the star cluster 
are always close to the adiabatic values if the stationary wind regime is possible. 

On the other hand, from the work of Chevalier & Clegg (1985), Canto et al. (2000), 
who considered the adiabatic solution for SSC winds, and from the recent semi-analytical 
radiative solution of Silich et al. (2004) we know that density and temperature inside of the 
star cluster volume remain almost constant (in the adiabatic case the ratio of the temperature 
at the star cluster surface to the central temperature is Tsc/Tc = 2/(7 + 1) = 0.75) and 
are completely defined by the star cluster parameters (Lsc, M S c and Rsc)- This is also the 
case in the radiative solution, although one has to take also the chemical composition of the 
ejected material into account. Thus the density and the temperature distributions in zone 
A may be approximated as 

T(r) w T sc , p{r) w a p p sc , (3) 

where the fiducial coefficient a p takes into consideration the deviation of density in zone 
A from the homogeneous distribution (the small deviation of temperature from the surface 
value leads to negligible changes of the X-ray emissivity). 

In zone B all variables rapidly approach their asymptotic values, p(r) ~ r~ 2 , T ~ r _4//3 . 
In the strongly radiative regime the temperature may drop faster and reach values ~ 10 4 K 
at smaller radii. However the contribution of this zone to the total X-ray luminosity is small 
and it arises mainly from a region close to the SSC surface (see section 4). 
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The expansion of the outer shock is supported by the thermal pressure of the thermalized 
wind (zone C in Figure 1). In the case of a homogeneous ISM and a constant rate of 
mechanical energy deposited by the central cluster, the time evolution of the outer shock 
radius, R ou t, the reverse shock radius, R in , and of the expansion velocity, V out , of the outer 
shell are given by (Koo & McKee, 1992) 

R^- o.88 1/5 1^, R m = om(^-) 3/W ^, v out = 3 -^, (4) 

where t is the bubble age and Pism is the density of the ISM. 

The shocked ISM in zone D remains adiabatic until the evolution time becomes larger 
than its characteristic cooling time scale, t coo i, (Mac Low & McCray, 1988), when 

t > t C ooi = ISZjgM x 10 4 n 7 sm L^yr, (5) 

where Zjsm is the ISM gas metallicity measured in the Solar units and hism is the ISM 
gas number density. When t exceeds t coo i, the gas in zone D becomes thermally unstable 
and collapses into a thin cold expanding outer shell. After this stage the density and the 
temperature distributions inside the bubble (zone C) are dominated by thermal evaporation 
of matter from the cold shell into the superbubble. If thermal evaporation is not blocked by 
magnetic fields, the density and the temperature distributions in zone C are (Weaver et al. 
1977; Mac Low & McCray, 1988) 

n(x) =n (l -x)- 2/5 ; T(x) = T (l - x) 2/ \ (6) 

where x = t/Rcd is the fractional radius, Rcd is the radius of the contact discontinuity, 
no and To are functions of the evolutionary time: uq = 1.0 x lO^gg 35 ^^ 5 ^ 22 / 35 cm -3 , 
T = 5.5 x 10 7 T3g 35 nj^t~ 6 / 35 K, and t here is measured in year units. 



4. The X-ray model 

Once the distributions of density and temperature in the various zones above described 
(zones A, B, C, including at earlier times also zone D) are known, one can calculate the 
X-ray luminosity from the whole remnant. This is given by 

rRout 

Tx = 4tt / r 2 n 2 A x (T,Z)dr, (7) 
Jo 

where n(r) is the atomic number density, R out marks the location of the outer shock and 
A X (Z,T) is the X-ray emissivity derived by Raymond & Smith (1977) in their hot-plasma 
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code (see Strickland & Stevens 2000). Clearly, regions in which the gas temperature is below 
the X-ray cut-off temperature (T cut rj 5 x 10 5 K), would not contribute to the global X-ray 
emission. 



4.1. The X-ray luminosity from the star cluster and its free wind 

The temperature and density distributions within zones A and B, obtained from the 
numerical models (see Silich et al. 2004), lead to the corresponding X-ray luminosity. The 
results of the calculations for a 10 6 M clusters of 1 pc and 10 pc radii are presented in Figure 
2. The figure shows the contribution, €f — [L x ,a+b—Lx(R)]/ Lx,a+b, arising from concentric 
shells with inner radius, R, and outer radius equal to the X-ray cutoff radius, to the total 
X-ray luminosity, Lx,a+b, from zones A and B. Figure 2 shows that the contribution from 
the free- wind region (zone B) is negligible. The value of e/ for shells with radius R > Rsc 
never exceeds ~ 20% of the cumulative emission from zones A and B. This is mainly due 
to the rapid decline in density (oc r~ 2 ) in the free wind region. A good estimate of L x ,a+b 
can be obtained if one approximates the density and temperature distributions within zone 
A by constant values (see equations 3). The evaluation of integral (7) then yields 

4a 2 p irA x p 2 sc R 3 sc 

^X,A+B - r-j , (P) 

where \i a = 14/11toh is the mean mass per atom. This simple analytic formula also accounts 
for the contribution from zone B, and is in excellent agreement with the results of our semi- 
analytical calculations (see Figure 3) if the fiducial coefficient, a p xs 2.0. One can also obtain 
psc from the equation (3) and re-write equation (8) in the form 

L x ,a+b = 3.8 x 10 34 $(T 5C , Z sc ) - L \ = 3.25 x 10 34 $(T SC , Z sc ) ^8 , (9) 

where Tk e v is the temperature of the X-ray plasma measured in keV units. Equation (9) 
is normalized to Ax = 3 x 10 -23 erg cm 3 s -1 and the normalization function, §(Tsc, Zsc), 
depends on the plasma temperature, T sc , and its metallicity, Zsc ( see Figure 4). It predicts 
a quadratic dependence of the X-ray luminosity (zones A and B) on the star cluster mass 
(Lsc scales linearly with the mass of the cluster) and is similar to the scaling relation from 
Stevens & Hartwell (2003; note a misprint in their equation 10 where a star cluster wind 
scaling parameter should be proportional to see also Oskinova, 2005). 

Equation (9) shows how the X-ray emission from the central zones A and B depends on 
the main star cluster parameters. The strong dependence on Va i00 comes from the fact that 
close to the star cluster surface p w ~ Msc/RscVa,oo ~ L S c/Vaoo- This last relationship 
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R/R sc 



Fig. 2. — The distribution of the radiated energy along the central zones A and B. Solid and 
dotted lines display the results of the calculations for quasi-adiabatic (solid line) and strongly 
radiative (dotted line) outflows. The star cluster mass and wind adiabatic terminal speed are 
10 6 M Q and Va j00 = 1000 km s _1 in both cases. The star cluster radii are Rgc — 10 pc (solid 
line) and R S c = 1 pc (dotted line). L x ,a+b = 3.2 x 10 39 erg s _1 and L x ,a+b = 2.6 x 10 38 
erg s" 1 in the case of 1 pc and 10 pc cluster, respectively. 
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Fig. 3. — The X-ray luminosity of SSCs and their winds (zones A+B) as a function of 
Rsc- The solid lines result from numerical calculations and are compared to the analytical 
formula (9) for star clusters with masses 1O 6 M and 1O 5 M (dotted lines). The adiabatic 
wind terminal speed is Va,oo — 1000 km s _1 in both cases, and the ejected gas metallicity is 
solar. 
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Fig. 4. — The normalization function $(T, Z) plotted as function of the X-ray plasma tem- 
perature measured in keV units. Different lines are marked with the assumed gas metallicity. 
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leads to equation (9) as the X-ray luminosity (along with radiative cooling) is in direct 
proportion to the square of the gas number density. This implies that the wind adiabatic 
terminal velocity (or the ratio of the energy to the mass deposition rates) is the key parameter 
that defines the X-ray emission from the central zones. Note that equation (9) cannot be 
applied to very compact and massive star clusters above the threshold line (see Silich et al. 
2004; Tenorio-Tagle et al. 2005), clusters for which the stationary wind solution is inhibited 
by radiative cooling. 



4.2. X-ray luminosity from the bubble 

4-2.1. Adiabatic stage. 

During the early stages, when t < t cooh the distributions of density and temperature in 
zones C and D are close to homogeneous and one can evaluate integral (7) analytically: 

4tt 233 3A x (T c ,Z c )L 2 sc (t - ^) 2 
L x ,c = ^-A X (T C , Z c )n 2 c (R 3 CD - R 3 in ) = 2 t/4 , ff3 psT = 

7r Pa V A,oo\ rL CD n in) 

9.5 x 1Q34 ( A * ) Affll- ^ erg 8 -l f (1Q) 



3xlO- 2 V V^ 00 [l-(R in /R CD ) 3 } 
L _4tt (pismV A x (T D ,Z ISM )Rl ut 

X ' D 3 UJ 1 " (RCD/Rout) 3 

4-6 x 10- (^z») «):>, I V'V <<■!,* . (11) 

where ric and Zc are the shocked wind (zone C) number density and metallicity, Tc and 
T D are plasma temperatures in zones C and D, and t 7 is the evolutionary time measured 
in units of 10 7 yr. Thus, at the adiabatic stage the luminosity from the outer shell highly 
exceeds that from the hot bubble interior and dominates the total X-ray output. This is due 
to the much smaller plasma density in zone C (pc/pd ~ M sc t/ pi SM R 3 out ~ V^JV^ 1). 
However, if the metallicity of the ISM is not well below the Solar value, this stage is so short 
(less than 10 5 yr; see equation 5), that we omit it in our further consideration. 



4-2.2. Hot bubble with a cold outer shell. 



At later times the swept up interstellar gas cools down to values well below the X-ray 
cut-off temperature and forms a dense outer shell separated from the hot bubble interior by 
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the contact discontinuity (see Figure 1). After this stage, the X-ray luminosity of zone C 
has been calculated by Chu et al. (1993). It is 

L X , C = 3.4 x l^Z c I(T)LftnTMtT*erg s~\ (12) 

where Zq is the metallicity of the X-ray plasma in zone C, I(r) is a dimensionless integral 
with a value close to 2. In the models of Chu et al. (1993) the main contribution to the 
bubble X-ray luminosity comes from the outer, densest layers of zone C, where the density 
and temperature are determined by thermal evaporation of the cold outer shell. For this 
reason the metallicity of the X-ray gas in this zone may be different from that in zones A 
and B and hardly ever exceeds a few times the solar value (see Silich et al. 2001). Hereafter 
we assume that the chemical composition of the plasma in zone C is solar, Zq = Z Q . On the 
other hand the temperature decreases towards the outer layers (see equation 6) and therefore 
the X-ray emission from zone C should be much softer than that from zones A and B. 

5. Integrated properties of the model predicted X-ray emission and 

comparison with observations 

5.1. Wind vs bubble luminosity 

We first used our analytic results (equation 9) to derive the expected X-ray luminosity 
from clusters of different masses and different radii. The calculated X-ray luminosities from 
zones A and B for a range of the star cluster wind terminal speeds (1000 km s _1 < Va i00 < 
1500 km s _1 ) are presented in Figure 5. The expected temperatures of the plasma are 
~ 1.1 x 10 7 K and ~ 2.5 x 10 7 K for V AoD = 1000 km s" 1 and V At<x = 1500 km s" 1 , respectively 
(see equation 1). The normalization functions are $ = 0.758 and $ = 6.52 for T ~ 1.1 x 10 7 K 
plasma and $ = 0.508 and $ = 3.16 for T ~ 2.5 x 10 7 K for a plasma with Z = Z Q and 
Z = 10Z Q , respectively (see Figure 4). Figure 5 shows that the X-ray luminosity from 
the cluster (zones A and B) is strongly dependent on the star cluster parameters and may 
easily vary within several orders of magnitude (the low luminosity limits move towards 
smaller values if Va jCO exceeds 1500 km s^ 1 ) even for a sample of clusters with a given 
mass. The crucial parameter, and unfortunately not a well known parameter, is the star 
cluster wind terminal speed, Va jCO = (2L SC / M S c) 1 ^ 2 ■ Perhaps the only reasonable way 
to have observational restrictions for this parameter for distant clusters is to measure the 
temperature of the X-ray plasma to then derive Va,oo from equation (1). This possibility 
may provide reasonable estimates for the most massive and compact clusters despite the 
presence of X-ray binaries (see below). 

The luminosities are also dependent on the metallicity of the X-ray plasma (compare 
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panels a and b in Figure 5). For more metallic plasmas the calculated luminosities are shifted 
to higher values and the difference between the low and the high velocity (temperature) 
outflows becomes even more significant as the enhanced metallicity affects less the emissivity 
of the higher temperature plasma (see Figure 4). 




Fig. 5. — The X-ray luminosity of SSCs and their winds (zones A+B) as a function of star 
cluster parameters. The shadow regions on the diagram represent (from the top to the bot- 
tom), the calculated X-ray luminosity for 1O 6 M , 10 5 M and 1O 4 M clusters, respectively. 
Va,oo = 1000 km s" 1 and Va,oo = 1500 km s" 1 were used to calculate the upper and lower 
luminosity limits for every star cluster. Panel a displays the results of the calculations for a 
plasma with solar metallicity, and panel b is the same for a plasma with Z=1OZ . 



We have also compared the X-ray emission from zones A and B with the X-ray emis- 
sion from the hot bubble interior (zone C) predicted by the Chu et al. (1993) standard 
model. Figure 6 shows the contributions from the star cluster (zones A and B) and from 
the expanding superbubble (zone C) to the total 0.3 keV - 8.0 keV emission, for clusters of 
different masses and radii evolving into an ISM with different densities. The X-ray emission 
from the high temperature zones A and B may dominate only in the case of very massive 
(> 1O 6 M ) and compact (few parsecs) star clusters evolving into a low density ISM (panel 
a in Figure 6). Otherwise the total luminosity is dominated by the hot superbubble interior. 
The contribution from the bubble plasma progressively increases with time and becomes 
dominant after a short while even for compact and massive clusters, if they evolve into a 
dense ISM (see panel b, Figure 6). The X-ray luminosity from star clusters decreases for 
larger clusters (see equation 9) shifting the luminosity from zones A and B well below the 
bubble luminosity even for massive clusters evolving into a low density ISM (see panel c in 
Figure 6). Clearly, this tendency becomes even stronger for less massive clusters (panels d, 
e, f in Figure 6). However, even in the case of the rather dense ISM (n IS M = 10 cm -3 ), the 
radius of the bubble exceeds 90 pc after 1 Myr, in the case of a 10 5 M cluster and 150 pc 
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Fig. 6. — X-ray luminosity generated by star clusters of different masses and radii located in 
different ISM environments, a) 1O 6 M 1 pc cluster in the low density ISM with riisu = 1 
cm -3 ; b) 10 6 M 1 pc cluster in the dense ISM with Uism = 10 cm -3 ; c) 10 6 M Q 10 pc cluster 
in the low density ISM with riisu — 1 cm" 3 ; d) 1O 5 M 1 pc cluster in the low density ISM 
with iijsm — 1 cm" 3 ; e) 1O 5 M 1 pc cluster in the dense ISM with Ujsm — 10 cm" 3 ; f) 
1O 5 M 10 pc cluster in the low density ISM with uism = 1 cm -3 ; Solid lines display the 
total X-ray luminosity. Dotted lines are contribution from zones A and B, and dashed lines 
are contribution from zone C, respectively. 
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in the case of a 1O 6 M cluster. This implies that compact and massive star clusters are 
to be detected as point-like X-ray sources embedded into extended regions of diffuse X-ray 
emission. 

One can also distinguish between zones A and B and zone C from their contributions 
to the hard and soft energy bands. Indeed, the temperature of the hot plasma ejected from 
the star cluster region (zone A) is defined by the ratio of the mechanical energy and mass 
deposition rates (or by the terminal speed of the star cluster wind) and, for reasonable 
values of Va j00 , falls in the range 1 keV - 10 keV (see equation 1). At the same time the 
temperature of the plasma which dominates the bubble X-ray emission is much lower due 
to dilution with the cold material evaporated from the outer shell. This suggests that the 
soft component to the X-ray emission is dominated by the bubble plasma whereas the hard 
component is defined by the hot thermalized ejecta (zone A and B). This conclusion is 
illustrated by Figure 7 which presents the results of the numerical calculations for a 10 5 M 
cluster of 1 pc radius evolving into a 1 cm~ 3 ISM (shown also in Figure 6d). Figure 7 clearly 
demonstrates that the star cluster ejecta dominates the hard X-ray emission even if the total 
luminosity is completely defined by the bubble plasma. Note that because of the assumed 
constant emissivity of the X-ray plasma, formula (12) overestimates the X-ray emission from 
the bubble interior by a factor ~ 2. 




Zone C 



Time (Myr) 



Time (Myr) 



Fig. 7. — Soft and hard components of the X-ray luminosity generated by the SSC. The 
contribution of different zones to the soft (0.3 keV -2.0 keV) and hard (2.0 keV - 8.0 keV) 
X-ray emission is displayed in panels a and b, respectively. The star cluster mass is 10 5 M , 
radius 1 pc, and the star cluster wind terminal speed is 1000 km s -3 . The density of the 
ISM is 1 cm" 3 . The metallicity was assumed to be Z=Z in all zones. 



From equation (9) one can derive the fraction, ex-, of the injected mechanical energy, 
L S c, that is transformed into the X-rays. The X-ray production efficiency, €a+b, for the 
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hard component associated with the star cluster plasma is 

e A+ B = = 3.8 x 10- 4 $(T 5C , Z sc ) L ™ . (13) 

L sc -Ksc.i^iooo 

For 10 6 M clusters with L 38 = 3 x 10 2 it hardly exceeds 10% even for most compact clusters 
and decreases progressively for star clusters with smaller masses. 

One can obtain the X-ray production efficiency, ec, for the soft bubble plasma from 
equation (12): 

e c = 3.4 x 10- 2 Z c /(r)L- 8 2/35 nS4 9/35 , (14) 

It is weakly dependent on the star cluster parameters, but is dependent on the density of 
the ISM and becomes larger as the evolutionary time grows, ec was discussed by Cerviho et 
al. (2002) who finally chose an arbitrary value of ec = 20%, and by Smith et al. (2005) who 
collected published data for well-studied HII regions and superbubbles and concluded that 
for young objects ec ~ 0.02%, while for older clusters ec falls into the range 0.2% - 7%. 



5.2. Comparison with nearby clusters 

A comparison of the model predictions with observed X-ray luminosities requires that 
parameters of the embedded star cluster (half-light radius, Rsc, mass, M$c and age, Tsc) will 
be obtained from the optical or infrared observations and then be compared with the observed 
X-ray luminosity, temperature (which is proportional to the star cluster wind terminal speed, 
Va,oo) and metallicity of the X-ray plasma. The full set of required parameters are available 
only for a restricted sample of nearby clusters. Below we confront our model (equation 
9) with several nearby young stellar clusters. In contrast with Stevens & Hartwell (2003), 
we will use for the comparison an observed temperature of the X-ray plasma which is an 
appropriate parameter for distant clusters, and will avoid to use the mean weighted terminal 
speed derived from the analysis of the individual massive stars, embedded members of the 
cluster. 



5.2.1. The sample of clusters 

A) NGC 3603. This cluster appears to be one of the densest and more massive cluster 
known in our Galaxy. Within a 1 pc radius NGC 3603 reveals a remarkable similarity in the 
central density and stellar density distribution of R136, the dense core of 30 Dor in the Large 
Magellanic Cloud (Moffat et al. 1994). However outside Rsc ~ 1 P c , the stellar density in 
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NGC 3603 decreases sharply whereas the stellar density in 30 Dor remains to decrease with 
a similar slope to larger radii. Malumuth & Heap (1994) measured the ionizing flux from 
30 Dor as a function of distance from the cluster center. They found that the number of 
ionizing photons from the inner 1 pc (4") region is N° fa 1.7 x 10 51 s _1 (see their figure 16). 
This corresponds to a 1.7 x 1O 4 M cluster if compared with Leitherer & Heckman (1995) 
evolutionary synthesis model and approximately two times smaller than the ~ 3 x 1O 4 M 
value obtained by Brandl et al. (1996) for the inner 20" region of 30 Dor. Because of the 
identical mass distribution in R136 and NGC 3603 inside the inner 1 pc radius, we adopt 
1.7 x 1O 4 M value as a total mass of the NGC 3603 cluster. Sung & Bessel (2004) used 
optical color-magnitude diagrams for massive members of NGC 3603 and concluded that 
the age of the cluster is t$c = 1 ± 1 Myr. The X-ray emission from the cluster has been 
studied by Moffat et al. (2002) who detected around 40 point-like X-ray sources definitely 
associated with the star cluster and a local diffuse emission most probably associated with 
the hot plasma inside the star cluster and with the free wind region. Moffat et al. (2002) 
concluded that approximately 20% (2 x 10 34 erg s _1 ) of the total 1 x 10 35 erg s _1 X-ray 
luminosity in the 0.5 — 10 keV band is associated with the local diffuse component of solar 
abundance and the best fitted temperature of kT core fa 3.1 keV at the core and kT fa 2.1 
keV (which is close to the temperature at the star cluster edge that is expected from the 
star cluster wind model) in the outer region. 

B) Arches cluster. The Arches cluster is another young (t S c = 2 ± 1 Myr; Figer et al. 
1999a) and very compact (Rsc ~ 0.2 pc; Figer et al. 1999b) Galactic cluster. The mass 
of the cluster has been estimated by Figer et al. (1999a) who extrapolated the observed 
number of massive stars down to the low cutoff mass limit. If a 1 M low mass cutoff is 
adopted, the star cluster mass would be Msc ~ 1.1 x 1O 4 M . The X-ray luminosity from 
the cluster was first detected by Yusef-Zadeh et al. (2002) who revealed within the volume 
occupied by the star cluster two point-like sources embedded into a more extended (down 
to a 1.2-1.8 pc radius) X-ray halo. Law & Yusef-Zadeh (2004) found that the observed X- 
ray emission may be equally well fitted with an overabundant, Z = (4 — 5)Z , one or two 
temperature thermal model, or with the one temperature plus power-law model with ~ 1/6 
of the flux attributed to the power-law component. The required temperature of the thermal 
plasma is around kT = 1.5 keV. The models give L x = (0.5 — 1) x 10 35 erg s _1 for the 0.5 
- 8.0 keV energy band. Overall, the point-like sources inside the cluster contribute ~ 60% 
to the total emission and the rest is distributed throughout the region with dimensions of 
approximately 3.6 x 2.4 pc. This leads to the X-ray luminosity from the cluster plasma, 
Lx = (2 — 4) x 10 34 erg s _1 , that is slightly larger than that from the earlier estimates 
of Yusef-Zadeh et al. (2002; L x = 1.6 x 10 34 erg s _1 ). However the interpretation of the 
observed X-ray emission is ambiguous. The observed spectrum may also be fitted by the 
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combined power law plus a 6.4 keV Gaussian component that models a fluorescent Fe K a 
emission from a nearby molecular cloud. If this is the case, the contribution of the fluorescent 
component may confuse the estimated X-ray emission from the cluster plasma. 

C) Quintuplet cluster. The Quintuplet cluster is the smallest and the oldest (t S c = 
4± 1 Myr, Figer et al. 1999b) one. The mass of the cluster is M sc = 8.8 x 10 3 M if a 1 M 
low mass cutoff is used (Figer et al., 1999a), and the radius is about of 1 pc (Figer et al. 
1999b). The X-ray emission from the cluster was detected by Chandra. After subtraction of 
the four point-like sources, the X-ray spectrum of the Quintuplet cluster is perfectly fitted 
by the thermal plasma model with the plasma temperature 2.4 keV at the star cluster center 
and a solar metallicity (Law & Yusef-Zadeh, 2004). The absorption-corrected halo emission 
in the 0.5 - 8 keV energy band is 1 x 10 34 erg s -1 (Law & Yusef-Zadeh, 2004). From the 
observed central temperature we derive a plasma temperature at the star cluster surface (see 
section 3), kT = 0.75 x 2.4 = 1.8 keV, and then use it as an input parameter for our model. 

5.2.2. The comparison with model predictions 

In order to compare the observed X-ray luminosities with that predicted by the star 
cluster wind model, we take an advantage of our equation (9) and combine it with input 
parameters obtained from the optical examination of the host clusters. To calculate the 
energy deposition rate, L sc , provided by multiple stellar winds and supernovae explosions 
inside the cluster, we used a starburst synthetic model of Leitherer & Heckman (1995) and 
scale it with the star cluster mass. For a standard, 1O 6 M , cluster we adopt a 10 40 erg s _1 
energy deposition rate before the first supernova explosion (t < 3 Myr), and approximately 
three times larger average value after the beginning of the supernovae era. The star clusters 
parameters, the observed and the calculated X-ray luminosities are collected in Table 1. 

The calculated X-ray luminosities are listed in the last column of Table 1. They agree 
Table 1: Star clusters parameters and X-ray luminosities 



Cluster 


M sc 


Rsc 


Age 


Lsc 


T &s 


Z 


Lx,obs 


Lx,mod 




(M ) 


(pc) 


(Myr) 


(erg s _1 ) 


keV 


(Z ) 


(erg s -1 ) 


(erg s" 1 ) 


NGC 3603 


1.7 x 10 4 


1 


1 ±1 


1.7 x 10 38 


2.1 ±0.7 


1 


2.0 x 10 34 


5.2 x 10 33 


Arches 


1.1 x 10 4 


0.2 


2±1 


1.1 x 10 38 


1.5 ±0.2 


4 


(2.0 - 4.0) x 10 34 


9.3 x 10 34 


Quintuplet 


8.8 x 10 3 


1 


4± 1 


4.4 x 10 38 


2.4 ±0.5 


1 


1.0 x 10 34 


1.2 x 10 34 
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with the observed values within a factor of 4 and do not show a systematic shift to the smaller 
of larger values relative to the observed luminosities. Clearly, the analysis of a larger sample 
of clusters that includes more massive objects is required to avoid uncertainties related with 
estimates of masses, ages and errors that result from the deviation of star distributions from 
the idealized initial mass function. For example, according to Leitherer & Heckman (1995) 
the energy input rate from a coeval star cluster changes sharply in a short time interval 
between 2 Myr and 4 Myr. A small mistake in the age estimate then may change the 
predicted X-ray luminosity within a factor of 10 or more. A smaller upper mass cutoff leads 
to the same uncertainty. For a ~ 10 4 M cluster the relative dispersion of the observed X-ray 
luminosity due to the discreteness of the stellar population, a x = [{Lx — Lx) 2 } 1 / 2 /Lx, may 
reach ax ~ 3 (see Figure 2 from Cerviho et al. 2002), and thus the Lx detected from a 
particular cluster may deviate from the mean value within a factor of ~ 3. It seems therefore 
that more massive and more evolved star clusters are better candidates to compare with our 
model predictions. 

Several more complications should be taken into consideration when comparing the 
results from equation (9) to the observed X-ray emission. In all young stellar clusters we have 
considered, the observed X-ray emission presents two components: a number of unresolved 
point-like sources, probably associated with the most luminous individual stars and local 
density enhancements resulting from the collisions of nearby stellar winds (see Ozernoy et 
al. 1997, and the results from Raga et al. 2001 and Rockefeller et al. 2004). The contribution 
from the point-like sources varies from ~ 20% in the Quintuplet cluster to ~ 80% in NGC 
3603 and cannot be distinguished from the star cluster diffuse component in distant clusters. 

After ~ 3 Myr the most massive stars leave the main sequence and explode as super- 
nova. This results in the formation of X-ray binary systems. The high mass X-ray binaries 
(HMXBs) contain a relativistic remnant of the supernova explosion (a neutron star or a black 
hole) which captures and accretes a wind from the secondary (O or B type), massive star. 
To transform a gravitational energy into the X-ray emission effectively, it is essential that 
an accretion disc is formed or a very strong magnetic field is present (see, for example, van 
Bever & Vanbeveren, 2000, and references therein). Thus, it takes 4-5 Myr for the HMXBs 
to become active, depending on the upper mass limit occurring in the cluster. The HMXBs 
phase is restricted by the life-time of the massive secondaries and is typically few times 10 7 
yr. In the low mass X-ray binaries (LMXB) the optical component is a low mass post-main- 
sequence star that fills its Roche lobe. In the instantaneous star formation LMXBs become 
active at later times and certainly dominate the X-ray emission from older clusters for which 
the SNe activity has terminated and the X-ray emission from zones A and B has vanished. 
LMXBs are often detected in globular clusters (see, for example, Maccarone et al. 2003, 
who found that in the NGC 4472, approximately 40% of the LMXBs are associated with 
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the globular clusters and approximately 4% of the globular clusters contain LMXBs). In the 
Milky Way around 10% of the LMXBs are associated with globular clusters which contain 
less than 1% of the stellar mass of the Galaxy (Liu et al., 2001). NGC 3603, Arches and 
Quituplet clusters are too young and do not contain active binaries. 

The population synthesis models for massive star clusters which take into consideration 
binary systems have been developed by a number of groups (see, for example, Iben et al. 
1995; Mas-Hesse & Cerviho, 1999; Dalton et al. 1995). Van Bever & Vanbeveren (2000) 
have used such a model to calculate the X-ray luminosity from the binary population of 
instantaneous starbursts with different fractions of binaries at the birth of the cluster. The 
HMXB component begins to contribute to the X-ray emission after « 4 Myr. The normalized 
per one solar mass X-ray emission does not depend on the initial mass of the cluster and 
soon reaches a value that varies between 10 32 erg s -1 Mg 1 and 10 33 erg s" 1 Mq 1 (see figures 
1 and 2 from their paper). 

In contrast to HMXBs, the normalized per unit solar mass X-ray luminosity from the 
thermalized star cluster plasma scales linearly with the star cluster mass: 

L XA+B /M SC = 2.9 x 10 33 $(T SC , Z S c)M 6 /Rsc,iT* eV , (15) 

where M 6 is the star cluster mass measured in units of 10 6 M Q . Figure 8 shows that for the 
most massive and compact young stellar clusters the X-ray luminosity of the thermalized star 
cluster plasma (zones A and B) may be comparable or even exceed that from the HMXBs. 

6. Conclusions 

We have developed a simple way of estimating the X-ray emission generated by super 
star clusters. The derived expression (equation 9) takes into consideration the intrinsic 
parameters of SSCs and their winds (Rsc, L$c, Va,ooj Zsc)- I n particular, we have found 
that the expected X-ray luminosity from SSC is highly dependent on the star cluster wind 
terminal speed (Lx ~ KToo)> a quantity, related to the temperature of the hot thermalized 
plasma within the SSC volume, and scales quadratically with the star cluster mass. The 
proposed relation seems to be in reasonable agreement with parameters of nearby clusters 
and their detected X-ray emissions. 

We have also compared the X-ray luminosity from the SSCs with the luminosity of 
the interstellar bubbles generated from the mechanical interaction of the star cluster winds 
with the ISM. We found that the soft component and the total X-ray emission are usually 
dominated by the superbubble plasma. The contribution from the SSC plasma may dominate 
only in the case of very massive (> 10 6 M ) and compact (few parsecs) star clusters evolving 
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Fig. 8. — X-ray luminosity of zones A and B per unit stellar mass as a function of the star 
cluster mass. The calculations have been performed for a 1 pc clusters with 1000 km s _1 
(Tsc = 0.95 keV) and 1500 km s^ 1 (T$c = 2.1 keV) terminal speeds indicated by the solid 
and dashed lines, respectively. The metallicity of the X-ray plasma is Solar. The shadow 
region on the diagram represents the X-ray luminosity range expected from the HMXB 
population (see van Bever & Vanbeveren, 2000). 
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into a low density ISM. However the hard (2.0 keV - 8.0 keV) component of the X-ray 
emission is usually dominated by the hot compact regions associated with the SSCs. This 
implies that compact and massive star clusters should be detected as hard compact X-ray 
sources embedded into extended regions of soft diffuse X-ray emission. 

On the other hand the comparison with the population synthesis models which take 
into consideration binary systems shows that the X-ray emission from the thermalized star 
cluster plasma normalized per unit stellar mass increases linearly with the mass of the cluster 
and may be comparable or even exceed that from the population of the HMXBs for the 
most compact and massive star clusters. Thus the thermalization of stellar winds and SNe 
ejecta, particularly in massive, young and compact super star clusters, may present an X-ray 
production that is comparable or even larger than that expected from the HMXB population. 
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